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Abstract 

We  find  an  optical  Raman  lattice  without  spin-orbit  coupling  showing  chiral  topological  orders  for 
cold  atoms.  Two  incident  plane-wave  lasers  are  applied  to  simultaneously  generate  a  double-well 
square  lattice  and  periodic  Raman  couplings,  the  latter  of  which  drive  the  nearest-neighbor  hopping 
and  create  a  staggered  flux  pattern  across  the  lattice.  Such  a  minimal  setup  can  yield  the  quantum 
anomalous  Hall  effect  with  a  large  gap-bandwidth  ratio  in  the  single  particle  regime,  while  in  the 
interacting  regime  it  achieves  the  Ji~J2-K spin  model,  with  the  nearest-neighboring  ( Jx )  and  next- 
nearest-neightboring  (/2)  exchange  coupling  coefficients,  and  the  three  three-spin  interacting 
parameter  ( K)  is  controllable.  We  show  that  the  Ji~J2-K spin  model  may  support  a  chiral  spin  liquid 
phase.  It  is  interesting  that  the  quantum  anomalous  Hall  state  can  be  detected  by  only  measuring  the 
Bloch  states  in  the  two  symmetric  momentum  points  of  the  first  Brillouin  zone.  Further,  we  also  show 
that  heating  in  the  present  optical  Raman  lattice  can  be  essentially  reduced  compared  with  the 
conventional  laser-assisted  tunneling  schemes.  This  suggests  that  the  predicted  topological  states  be 
reachable  with  the  current  experimental  capability. 


1.  Introduction 

Generation  of  synthetic  gauge  fields  for  cold  atoms  opens  a  new  direction  in  the  study  of  exotic  topological  states 
beyond  natural  conditions.  Two  different  scientific  paths  have  been  followed  in  the  experiment  to  create 
synthetic  gauge  fields  via  optical  means.  One  is  to  adopt  Raman  couplings  between  different  internal  hyperfine 
levels  (atomic  spins)  [1-6],  which  has  recently  been  used  in  experiments  to  generate  synthetic  spin-orbit  (SO) 
coupling  for  cold  atoms  [7-10] .  Another  is  to  adopt  laser-assisted  hopping  between  neighboring  lattice  sites 
without  spin  flip,  which  can  generate  U(l)  fluxes  by  imprinting  the  phases  of  Raman  lasers  into  the  hopping 
matrix  elements  [11-15].  Compared  with  the  technique  using  spin- flip  Raman  couplings,  the  latter  strategy  can 
be  achieved  with  far- detuned  lasers,  and  therefore  can  avoid  the  spontaneous  decay  of  excited  states. 

Realization  of  a  gapped  insulating  topological  state  typically  necessitates  an  optical  lattice  and  synthetic 
gauge  fields  which  satisfy  proper  conditions  [12-28].  In  the  conventional  techniques,  the  optical  lattice  and 
gauge  fields  are  generated  through  different  atom-laser  couplings.  In  such  cases  the  topological  regimes  are 
achieved  with  careful  manipulations  of  parameters,  which  might  be  challenging  for  the  experimental 
observation.  Recently,  it  was  proposed  that  creations  of  the  optical  lattice  and  SO  couplings  can  be  integrated 
through  the  same  standing  wave  lasers,  and  this  new  technique  can  have  explicit  advantages  in  realizing 
topological  phases  with  minimal  setups  and  without  complicated  manipulations  |  1, 29].  Nevertheless, 
generating  SO  couplings  requires  near- resonant  light  which  heats  up  the  system  by  spontaneous  emission 
[7-10].  A  possible  resolution  of  this  difficulty  is  to  consider  lanthanide  atoms  which  can  have  less  heating  due 
to  large  fine  structure  splitting  and  narrow  natural  linewidth  in  the  excited  levels  [30]. 


©  20 1 6  IOP  Publishing  Ltd  and  Deutsche  Physikaiische  Geseiischaft 


IOP  Publishing 


New  J.  Phys.  18  (2016)  035004 


X-JLiu  etal 


Figure  1 .  (a)  The  laser  beam  incident  from  x  direction  has  nonzero  polarization  components  along)/  and  z  axes,  with  the  z  polarization 
field  having  two  components  of  wave-vectors  A)  and  2/q.  With  the  reflections  of  the  three  mirrors  Mj  (j  =  1,  2,  3),  a  double-well 
lattice  is  generated  and  can  be  controlled  by  the  1  / 4-wave  plate  and  electro-optic  modulator  (EOM);  (b)  the  created  double-well  lattice 
has  an  energy  off  set  between  A  and  B  sublattices;  (c)  an  additional  linearly-polarized  running  laser  beam,  with  polarization 
components  in  x  and y  directions,  is  applied  along  z  direction.  This  beam,  together  with  the  laser  components  used  to  create  square 
lattice,  can  induce  spatially  periodic  Raman  couplings  (Vrx,  Vr^)  between  A  and  B  sites,  as  illustrated  in  (d).  The  Raman  transitions 
also  generate  a  staggered  flux  pattern  simultaneously  for  the  nearest-neighbor  hopping. 


In  this  paper,  we  introduce  the  model  of  optical  Raman  lattice  without  SO  coupling  to  observe  chiral 
topological  phases  for  cold  atoms.  The  setup  includes  a  double-well  square  lattice  and  periodic  Raman  couplings 
generated  simultaneously  through  two  incident  plane- wave  beams.  We  show  that  this  scheme  can  naturally 
realize  chiral  topological  phases  without  fine  tunings,  and  may  have  advantages  in  the  experimental  observation 
including  the  minimized  heating  and  full  controllability  in  parameters. 

The  manuscript  is  organized  in  the  following  way.  In  section  2  we  present  the  model  realization,  and  discuss 
the  properties  of  the  model;  in  section  3  we  give  the  tight-binding  model,  which  shows  the  quantum  anomalous 
Hall  (QAH)  effect;  in  section  4  we  show  that  the  QAH  states  can  be  detected  by  measuring  only  Bloch  states  at 
two  symmetric  momenta  of  the  first  Brillouin  zone  (FBZ).  The  heating  effect  is  discussed  in  section  5,  and  the 
chiral  spin  liquid  phase  is  studied  in  section  6.  Finally,  we  present  the  conclusions  in  section  7. 

2.  Model 

In  this  section  we  present  the  details  of  the  model  realization.  The  system  includes  a  double-well  square  optical 
lattice  generated  by  an  incident  plane-wave  light  with  both  in-plane  and  out-of  plane  polarized  components, 
and  the  two-photon  Raman  potentials  which  induce  hopping  between  nearest-neighboring  sites  and  create  a 
staggered  flux  pattern  across  the  double  well  lattice.  The  lattice  and  Raman  potential  profiles  exhibit  the  relative 
antisymmetry  in  the  spatial  inversions. 

We  first  introduce  the  generation  of  a  2D  double  well  square  lattice  depicted  in  figures  1  (a)  and  (b),  with  an 
onsite  energy  difference  A  between  A  and  B  sites.  Based  on  the  experiments  by  NIST  group  [31],  this  lattice 
configuration  can  be  realized  via  an  incident  plane-wave  laser  beam  which  has  both  nonzero  in-plane  and  out- 
of-plane  linearly  polarized  components,  and  with  the  assistance  of  three  mirrors  [see  figure  1(a)].  The  total 
electric  field  of  the  incident  laser  beam  can  be  described  as 

E(x)  =  E0(cos  ay  +  sina^e1^1*-^)  +  £0zel2(-klX~ut\  (1) 

where  the  polarization  angle  a  determines  the  magnitudes  of  in-plane  and  out-of-plane  polarization 
components  with  wave- vector  /q,  and  the  out-of-plane  polarized  field  has  another  component  with  its  wave 
vector  (2ki)  being  twice  of  that  of  the  former  ones.  Note  that  all  the  components  of  the  incident  beam  can  be 
generated  from  a  single  laser  source  through  an  optical  frequency  doubler,  with  which  one  can  control  the  ratio 
of  the  field  strengths  \E0/E0\  in  the  experiment  [32].  With  the  reflection  by  mirrors  the  y -polarization 
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component  of  the  incident  laser  beam  changes  to  be  x  -polarization  component  when  the  laser  beam  propagates 
along  d=y  direction,  while  the  z- polarization  components  do  not  change  for  the  entire  optical  paths.  The  in¬ 
plane  polarized  components  ( x  and  y  polarization  components)  generate  the  standing  wave  as 

E xy  —  2 E0  cos  a  [  cos(kix)y  +  cos(kiy)x]  e~luJt.  (2) 

Here  we  have  neglected  all  irrelevant  constant  phase  factors  which  have  no  effect  on  our  results  (see  appendix  for 
details).  On  the  other  hand,  the  out-of-plane  light  components  can  interfere  and  generate  the  standing  wave  as 


Ezz  =  4Eq  sin [k\ (x  +  y)]sm[ki(x  -  y)]e 


— i2  uitz 


+ 


4 E0  since  sin +  /)  jcosj^y 


(x-y) 


(3) 


Furthermore,  we  consider  here  the  blue-detuned  optical  dipole  transitions.  The  total  lattice  potential  reads 

Vsq(x,  y)  =  y0[cos2(/cix)  +  cos2(/ci/)] 

+  Vq  sin2 [ki (x  +  y)]sin2[ki(x  -  y)] 

+  A  sin2|~-y(;c  +  y)  jcos2  \~(x  ~~  y)  • 


(4) 


The  amplitudes  are  taken  that  A  <  Vq,  and  the  above  potential  describes  a  double- well  square  lattice  illustrated 
in  figure  1  (b),  with  the  staggered  onsite  energy  offset  A  between  A  and  B  sites  well  controlled  by  the  polarization 
angle  a.  When  A  is  large  compared  with  the  bare  hopping  couplings  between  neighboring  s-orbitals  of  the  A  and 
B  sites,  the  effective  tunneling  between  them  is  suppressed,  while  the  diagonal  AA/BB  hoppings  (denoted  by 
t'A/B)  are  allowed  along  dashed  lines  in  figure  1  (b).  The  second  term  in  equation  (4)  reduces  the  difference  in 
height  of  the  barriers  along  the  AB-bond  and  the  diagonal  (AA/BB)  directions.  Thus  it  can  enhance  tfA/B  relative 
to  the  hopping  coupling  between  A  and  B  sites,  providing  vast  tunability  in  parameters. 

The  tunneling  between  neighboring  A  and  B  sites  (denoted  by  tjj)  can  be  restored  by  two-photon  Raman 
couplings.  A  key  ingredient  of  the  present  scheme  is  that  the  in-plane  blue-detuned  laser  beam  which  generates 
the  square  lattice  also  takes  part  in  the  generation  of  Raman  couplings.  For  this  we  apply  an  additional  plane- 
wave  laser  beam  with  frequency  a;  —  6u  (Suj  ~  A),  propagating  along  the  perpendicular  z  direction  and  having 
linear  polarization  components  along  x  and  y  axes  (figure  1  (c)).  This  beam  is  described  by 

£  xy(z)  =  EiCe^x  +  e^rp)ezH  (5) 

where  (j)xjy  is  the  initial  phase  of  the  x/y- axis  polarization  component.  With  the  assistance  of  both  Exy  and  Exyy 
two  independent  Raman  couplings  are  induced  by  the  x-  and  y -polarization  components,  respectively 
(figure  1(d)).  In  particular,  the  x  (j>)- components  of  the  lights  Exy  and  Exy  generate  the  Raman  potential 
Vrx  (VRy)  which  takes  the  form 

VkxVty)  =  VR  cost^OOle^+^w  +  C.C.,  (6) 

where  the  amplitude  VR  oc  E0Ei  cos  a.  We  shall  see  below  that  a  finite  magnitude  of  <fix  —  cj)y ,  controllable  in 
experiment,  gives  rise  to  a  nonzero  staggered  flux  pattern  for  the  square  lattice,  as  illustrated  in  figure  1  (c). 

From  equations  (4)-(6)  we  can  see  that  the  zeros  of  VRx,Ry  are  located  at  the  lattice-site  centers,  which  implies 
that  the  Raman  potentials  are  parity  odd  relative  to  each  lattice- site  center  (figure  1  (c)).  With  this  key  property 
the  present  blue-detuned  optical  Raman  lattice  can  naturally  realize  topological  states  and  exhibit  essential 
advantages  in  minimize  heating  effects  for  experimental  studies.  The  symmetry  properties  of  s-orbitals  and 
VRXyRy  lead  to  two  important  consequences.  First,  the  Raman  potential  VRx  (VRy)  only  induces  the  nearest- 
neighbor  hopping  along  x  (y)  direction.  The  hopping  along  x/y  axis  is  associated  with  a  phase  (j)y/x  (—  <fiy/x)  if  the 
hopping  is  toward  (away  from)  B  sites.  In  experiment,  one  can  set  that  </>x  —  <fiy  —  20o,  which  is  equivalent  to  put 
that  <j)x  =  —  (f)y  =  <fi0.  Then  the  hopping  along  the  directions  depicted  by  arrows  in  figure  1  (c)  acquires  a  phase 
0O,  resulting  in  a  staggered  flux  pattern  with  the  flux  <F  =  4 cf)0.  Secondly,  the  hopping  from  one  site  to  its 
leftward  (upward)  neighboring  site  has  an  additional  minus  sign  relative  to  the  hopping  to  its  rightward 
(downward)  neighboring  site.  It  is  important  that  all  these  interesting  properties  are  obtained  automatically  by 
using  the  two  incident  beams  without  complicated  fine  tunings. 


3.  Quantum  anomalous  Hall  effect  with  large  gap-bandwidth  ratio 

Now  we  give  the  tight-binding  Hamiltonian.  Based  on  the  previous  analysis  we  can  obtain  the  s-band  tight- 
binding  Hamiltonian  in  the  following  form 
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Figure  2.  (a)  The  band  gap  £gap  versus  Raman  potential  amplitude  VR  (in  units  of  the  recoil  energy  ER)  and  the  hopping  phase  c/)0;  (b) 
the  ratio  f0/ t\  with  t\  =  0.5  (t'A  +  t'B)  and  (c)  the  ratio  £gap/£width  (blue  curve)  as  functions  of  VR.  The  phase  r/>0  =  7t/4  (n/8)  for  solid 
(dashed)  curves;  (d)  the  bulk  spectrum  with  VR  =  0.65ER  and  c/)0  =  7r/4.  Other  parameters  are  Vo /ER  =  4,  Vo/ER  =  l,and 
6uj  =  A  —  0.6Er.  For  87Rb  atoms  using  A  =  532  nm  yields  ER  «  27T  x  8  kHz,  and  t\  ~2i r  x  27  Hz.  The  Chern  number  Q  =  +1. 


H  =  -  cos  4>0  + wy  sin  (f)0) cl  7cAj 

(17) 

-  X  X  +  m*X(«M  - 

((i,J))H=A,B  ^  1 


where  ny  ^  =  J  fix j,  with  cj  j  and  c^j  the  creation  and  annihilation  operator,  and  the  Zeeman  term 

mz  =  ( A  —  6cj) /2.  The  notations  (z ,  j )  and  ( (i ,  j ) )  represent  the  summations  over  the  nearest-  and  next- 
nearest- neighboring  sites.  The  hopping  phase  (j)0  determines  the  real  (proportional  to  cos  </>0)  and  imaginary 
(proportional  to  sin  0O)  parts  of  the  nearest  hopping  coefficients,  with  the  staggered  factor  vj-  =  1  (  —  1)  for 
hopping  along  (opposite  to)  the  marked  direction  in  figure  1(c).  From  the  periodic  profiles  of  the  Raman 
potentials,  the  nearest-neighbor  hopping  coefficients  satisfy  tjj±lx  =  =t(— 1  )lxt0,  t^j± ly  =  =p(—  l)lxt^  with 

to  =  VR  J d2r^p  (r)  sin  (hx) V’^’s0)  (r) ,  (8) 

and  the  diagonal  hopping  coefficient 


t'r  =  J>  <;,< 


'(r) 


2m 


+  Vsq(r) 


^(f). 


(9) 


with  ^S(r)  ip  =  A,  B)  the  s-orbital  wave  function  at  the  j  -th  site.  We  note  that  the  staggered  sign  factor  ( —  l)1* 
is  due  to  the  periodicity  of  Raman  potentials  and  staggered  position  distribution  of  A  and  B  sites,  and  can  be 
absorbed  by  redefining  the  annihilation  operator  of  B  sites  to  be  cB  j  =  emxi/acBj,  with  a  the  lattice  constant.  In 
terms  of  the  new  basis,  the  diagonal  hopping  coefficient  for  the  B  sites  reverses  sign  tB  — ►  —tB.  The  tight-binding 
model  can  now  be  obtained  directly  and  in  k  space  the  Bloch  Hamiltonian  reads 
7Y(k)  =  —  2 10  cos  </>0(sin  kxa  +  sin  ky  a)  ax  —  2 10  sin  0O  (sin  kxa  —  sin  kya)ay 

+  [mz  —  2  (t'A  +  tB)  cos  kx  a  cos  ky  a\  az ,  with  the  Zeeman  term  mz  =  ( A  —  6u) / 2  and  crx^yZ  the  Pauli  matrices. 
It  is  interesting  that  without  Zeeman  and  diagonal  hopping  terms,  i.e.  if  mz  =  tfA  B  =  0,  the  above  Hamiltonian 
would  describe  massless  Dirac  fermions  with  two  independent  Dirac  points  at  A[  =  (0,  0)  and  A2  =  (0,  7 r). 
The  bulk  is  gapped  when  mz  ^±2  (tfA -\-  tB)andcj)0  ^  n7r/2  with  n  being  integer.  The  QAH  phases  [17, 33, 34] 
are  obtained  for  | mz \  <  l(tA  +  tB),  with  the  first  Chern  number  Q  =  sgn(</>0)  (0  <  |0O|  <  7r/2),  and  trivial 
regime  results  for  |  mz \  >2  {tA  +  tB).  Figure  2  provides  the  numerical  estimate  with  V0/ER  =  4,  V0/ER  =  1, 

A  =  0.6F^,  and  the  recoil  energy  ER  «  2tt  x  8  kHz  using  A  =  27r/fci  =  532  nm  for  87Rb  atoms,  which  gives 
thdLttA  B  ~2ttx27  Hz.  By  setting  0o  =  7r/4andmz  =  0,  the  bulk  gap  £gap  =  4  (tA+  tB)  ^  2i r  x  0.21kHz 
when  t0  >  tA  +  t^forV^  >  0. 7 1 ER  (figure  2 (a)).  Figures  2  (b)-(d)  show  a  large  ratio  (~4. 9)  between  the  band 
gap  and  bandwidth  £width  in  the  range  from  t0  =  0.7  (tA  +  tB )  at  Vr  ^  0.5  ITr  to  t0  =  tA  +  tB  at  Vr  ~  0.7Tr.  It 
is  noteworthy  that  a  large  gap-bandwidth  ratio  can  enable  the  study  of  correlated  topological  states  like  the 
fractional  QAH  effect  [35]  in  the  interacting  regime. 


4.  Detection  with  minimal  measurements 

Detection  of  the  QAH  insulating  phase  can  be  carried  out  with  several  different  measurement  strategies  in  the 
edge  [17,19]  and  bulk  [36-40] .  In  particular,  it  was  proposed  recently  that  the  topology  of  a  QAH  insulator  can 
be  determined  by  measuring  Bloch  eigenstates  at  only  two  or  four  highly  symmetric  points  of  the  first  Brillouin 
zone  [4 1].  It  was  shown  that  this  approach  is  valid  for  QAH  insulators  which  satisfy  the  inversion  symmetry 
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defined  by  P  =  P  0  R2d>  where  R2D  is  a  2D  spatial  inversion  operator  transforming  the  Bravais  lattice  vector 
R  — ►  —  R  and  P  is  a  parity  operator  acting  on  the  (pseudo)spin  space  [41].  In  the  present  lattice  system  the  unit 
cell  is  doubled  relative  to  the  original  square  lattice,  and  from  the  Hamiltonian  TC  (k)  one  can  check  that  no 
parity  symmetry  can  be  satisfied.  Nevertheless,  we  show  below  that  this  minimal  measurement  method  can  be 
still  applied  to  the  present  system  with  a  highly  nontrivial  generalization. 

In  the  physical  Hamiltonian  Tt  (k)  the  Pauli  matrices  <Jx,y,z  operate  on  the  sublattice  space.  To  complete  our 
proof,  we  construct  an  artificial  Hamiltonian  which  is  formally  equivalent  to  H  (k):  H(k)  =  dx  (k)  bx 
+dy  (k)  by  +  dz  (k)  bz .  The  only  difference  is  that  in  the  new  Hamiltonian  we  assume  that  bXyytZ  act  on  a  spin 
space  which  is  independent  of  the  position  space.  In  this  way,  we  know  that  the  new  Hamiltonian  7Y(k)  is 
invariant  under  the  following  2D  inversion  transformation  on  both  the  position  and  spin  space 

P  =  bz®  R2 d,  (10) 

where  the  transformation  R2D  sends  (kx,  ky)  — >  (—kx,  —ky).  Therefore,  at  the  four  symmetric  points 
(A  i)  =  {(0,  0),  (0,  7 r),  (7 r,  0),  (7 r,  7r) }  the  Bloch  states  are  also  parity  eigenstates  with 
P|w±(A*))  =  £fz)\u±(Ai)),  and  =  +1  or  —  1.  The  topology  of  the  artificial  insulating  system  can  be 
determined  by  the  following  invariant  [4 1  ] 

(-1)"  =  n  £<-)(A;) 

i 

4 

=  ]^[  sgn[mz  —  2(t'A  +  £B)cosAz,acosAz,a] 

i 

=  [(-ir]2  =  1.  (id 

In  the  third  line  of  the  above  equation  we  have  defined  that 

2 

(-D"  =  n  c<_)(Ai)  =  sgn[wz  -  4  (t’A  +  trB )2].  (12) 


Therefore  the  invariant  for  the  constructed  system  u  =  0.  This  indicates  that  the  Chern  number  for  the 
Hamiltonian  7Y(k)  should  always  be  even  Q  =  IN  [41].  This  result  is  easy  to  understand.  As  pointed  out 
previously,  for  the  original  physical  system  the  unit  cell  is  doubled.  Accordingly,  the  FBZ  of  the  original  square 
lattice,  denoted  by  DFBZ  is  only  half  of  the  FBZ  DFBZ  for  the  constructed  artificial  system.  The  two  momenta 
( kx ,  ky)  and  ( kx  +  7r,  ky  +  7r)  correspond  to  the  same  point  in  QFBZ  (note  that  k  is  not  the  sublattice 
momentum,  but  the  momentum  of  the  original  lattice  system  which  includes  both  A  and  B  sublattices). 
Therefore,  the  Chern  number  for  the  Bloch  Hamiltonian  H(k)  reads 

Q  =  —  f  dkxdkyn  •  (dkxn  x  d\  n) 

47T  dkeOpBz 

=  —  f  dkx dkyn  ■  ( dkxn  x  dk  n) 

Ait  dkeoFBZ 

+  —  f  dkxdkyn  ■  (, dkxn  x  dk  n) 

Ait  ^keQFBz+(7r^) 

=  A  f  dkxdkyn  ■  ( dkn  x  dk  n) 

27 T  dkGOpgZ 

=  2Q.  (13) 


From  the  number  v  one  cannot  tell  the  difference  of  a  topological  phase  from  a  trivial  phase.  In  the  next  step,  we 
shall  show  that  the  topology  of  the  artificial  system  can  also  be  determined  by  the  invariant  v  which  is  defined  in 
equation  (12)  with  the  parity  eigenvalues  at  Ai  =  {(0,  0)and  A2  =  (0,  7r),  half  of  the  four  parity- symmetric 
points  in  DFBZ.  The  magnitudes  v  =  0  and  +  1  correspond  to  the  topologically  trivial  and  nontrivial  states, 
respectively.  Then,  together  with  the  above  relation,  we  can  further  use  this  invariant  to  characterize  the  topology 
of  the  original  physical  system. 

The  proof  is  straightforward  and  is  valid  for  any  two-band  system  satisfying  the  following  two  conditions. 
First,  the  quantum  anomalous  Hall  phases  are  characterized  by  low  Chern  numbers.  In  particular,  for  the 
artificial  system  it  is  Q  =  {0,  ±2 }  and  for  the  original  physical  system  Q  =  {0,  d=  1 }.  Second,  the  system  can  be 
adiabatically  connected  to  the  one  obtained  under  a  four- fold  C4  rotational  transformation  on  such  system.  In 
other  words,  the  topology  is  not  changed  under  the  C4  transformation  in  position  and  (pseudo)spin  space 

MuHMp  ~  n,  M4  =  (14) 

where  R4  {it/ 2)  is  the  7r/2 -rotation  on  the  position  space,  transforming  the  Bloch  momentum 
{ kx ,  ky)  — >  ( ky ,  —  kx).  It  is  easy  to  see  that  the  inversion  symmetry  in  equation  ( 10)  is  given  by  P  =  A4j.  By  a 
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direct  check  one  can  verify  that  the  constructed  system  in  our  consideration  belongs  to  the  class  of  Hamiltonians 
satisfying  the  above  conditions.  What  we  need  to  prove  is  that  the  transition  between  a  trivial  phase  and  a 
topological  phase  must  be  associated  with  the  change  from  v  =  0  to  z/  =  + 1.  Let  the  system  be  initially  a  trivial 
insulator.  To  have  topological  phase  transition,  the  bulk  gap  should  close  and  reopen  at  some  momentum 
points.  Around  such  momenta  the  bulk  can  be  described  by  massive  Dirac  Hamiltonians,  with  the  Dirac  masses 
changing  signs  during  the  transition.  We  denote  one  of  the  Dirac  momentum  as  km  =  (fcb  k2).  Then,  from  the 
C4  symmetry  we  know  that  there  are  four-fold  of  such  Dirac  points  k Dj  (j  =  1 ,  2,  3,  4).  Moreover,  from  the 
relation  between  the  artificial  and  original  physical  systems,  we  have  that  the  momentum  k Dj  +  (7 r,  7 r)  is  also  a 
Dirac  point.  With  these  results  in  mind,  we  get  that  when  a  topological  phase  transition  occurs,  the  Dirac  masses 
simultaneously  reverse  signs  at  following  momenta  (not  necessarily  independent) 

km  =  (ki>  k2),  km  =  (k2,  —  ki),  km  =  (— kb  —  k2), 

kD4  =  (~k2,  ki),  kDj+4  =  k Dj  +  (7 r,  7 r),  j  =  1,...,4.  (15) 

It  is  easy  to  know  that  there  must  be  even  number  (denoted  as  2 N)  of  Dirac  points  in  the  above  formula  which 
are  independent.  On  the  other  hand,  from  the  symmetry  we  know  that  all  these  Dirac  points  contribute  the  same 
Chern  number  to  the  whole  bulk  invariant.  Before  and  after  the  phase  transition  the  Chern  number  changes  by 

f-4,  final  Cl,  initial  —  2N .  (Id) 

We  have  three  different  cases.  First,  if  km  is  an  inversion  symmetric  point,  e.g.  =  Ai  =  (0,  0),  the 
equation  (15)  includes  only  two  independent  points,  km  and  km  +  (77,  7r).  Second,  for  the  case  with 
km  =  (tt/2,  7t/2),  the  equation  (15)  includes  four  independent  Dirac  points,  i.e.  {kD;}  =  {(±7t/2,  d=7r/2)}. 
Finally,  for  the  remaining  cases  all  the  8  momenta  in  equation  (15)  are  independent.  This  implies  that  for  the 
later  two  situations  the  Chern  number  changes  by  4  and  8,  respectively,  while  in  the  first  case  the  Chern  number 
changes  by  2.  Therefore,  for  a  system  with  low  Chern  number,  only  the  first  situation  can  happen,  namely,  the 
bulk  gap  must  close  and  reopen  at  two  inversion  symmetric  momenta  A\  and  A3  (or  A2  and  A4).  Note  that  the 
Dirac  masses  at  these  points  are  equivalent  to  the  parity  eigenvalues,  so  the  topological  phase  transition  must  be 
associated  with  the  sign  change  of  corresponding  parity  eigenvalues,  leading  to  the  change  of  the  invariant  v. 
Furthermore,  it  is  easy  to  verify  that  the  trivial  phase  with  Q  correspond  to  v  =  0,  and  then  the  topological 
phases  with  Q  =  d=2  are  given  by  v  =  +1.  Together  with  the  relation  (13)  we  conclude  that  the  invariant  v 
classifies  the  topology  of  the  original  physical  system.  This  completes  our  proof. 

Since  the  parity  operator  is  crz,  the  parity  eigenvalues  are  the  pseudospin  eigenstates.  To  measure  the  parity 
eigenvalues  one  can  measure  the  pseudospin  polarization,  i.e.  the  population  difference  of  atoms  between  A  and 
B  sublattices,  which  can  be  measured  with  in  situ  imaging.  The  pseudospin  polarization  is  defined  by 
Ps  =  (Na  ~  Nb)/(Na  +  Ng),  where  NA/B  represents  the  number  of  atoms  in  the  A/B  sublattice.  It  follows  that 

2 

(-lr  =  n  sgn[p5(A®)].  (17) 

i=  1 

The  topological  phase  corresponds  to  v  —  —  lfor|m2|  <  2{t'A  +  tg),  and  the  trivial  phase  corresponds  to  v  =  0 
for\mz\  >  2(t'A  +  tg).  The  detection  can  be  carried  out  with  a  pseudospin- resolved  Bloch  oscillation  [41].  With 
an  external  force  applied  along  the  x  direction,  the  momentum  of  an  initial  atomic  cloud  evolves  along  the 
direction  from  Ai  =  (0,  0)to  A2  =  (0,  7r)  (figure  3).  For  the  topological  phase,  the  pseudospin  polarization 
ps  (t)  reverses  sign  from  sgn  [ps  (0)  ]  =  —  1  to  sgn  [ps  (r )  ]  =  + 1  at  half  Bloch  time  r  =  Tg/2  and  returns  to 
sgn[ps(r)]  =  —  latr  =  (figure  3(a)  and  (b)).  On  the  other  hand,  in  the  trivial  regime  \mz\  >  2  (tA  +  tg),  the 
sign  of  the  polarization  keeps  unchanged  during  the  Bloch  oscillation  (figure  3(c)  and  (d)).  Only  qualitative 
measurements  at  the  two  symmetric  momenta  are  needed  for  the  experimental  detection. 


5.  Heating 

The  experimental  feasibility  of  observing  topological  states,  especially  the  correlated  topological  states,  crucially 
depends  on  heating  effects  in  the  realization.  Note  that  in  the  present  scheme  all  applied  lights  are  far-detuned. 
Thus  the  heating  due  to  spontaneous  decay  from  excited  atomic  levels  is  negligible.  On  the  other  hand,  due  to  the 
energy  difference  between  the  Raman  photons,  the  two-photon  Raman  processes  have  two  main  consequences. 
The  first  one  is  that  it  induces  hopping  between  neighboring  sites  by  compensating  the  energy  difference 
between  A  and  B  sites,  leading  to  the  time-independent  effective  tight-binding  Hamiltonian  (7)  and  the 
topological  band  as  discussed  in  the  previous  sections.  This  process  includes  both  the  energy  transfer  from 
Raman  photons  to  atoms  (hopping  from  B  to  A  sites)  and  energy  transfer  from  atoms  to  Raman  photons 
(hopping  from  A  to  B  sites).  Thus  in  average  it  does  not  have  net  energy  transfer  between  Raman  photons  and  the 
atoms.  Another  main  consequence  is  that  if  the  two-photon  Raman  process  drives  an  onsite  transition  rather 
than  the  hopping,  it  can  transfer  the  energy  from  Raman  photons  to  atoms.  This  process  leads  to  the  increase  of 
the  total  mean  mechanical  energy,  i.e.  the  expectation  value  of  the  kinetic  and  potential  energies  [42],  of  an  atom 
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Figure  3.  The  time  evolution  of  numbers  of  atoms  populated  in  A  sites  (NA,  in  red  curves)  and  B  sites  (NB,  in  black  curves),  and  the 
polarization  ps  (r)  (blue  curves)  in  the  Bloch  oscillation,  with  N0  —  NA  +  NB.  (a)  and  (b)  are  topological  regimes  with  mz  —  0  and 
mz  =  —l.5(t'A  +  tB),  respectively;  (c)  and  (d)  are  trivial  regimes  with  mz  =  —2.5(tA  +  tB)  and  mz  —  —3  (tA  +  tB),  respectively.  The 
hopping  phase  is  taken  as  cf>0  =  7t/4. 


trapped  in  the  lattice  potential.  Note  that  the  laser- assisted  tunneling  scheme  without  spin  flip  does  not  suffer 
from  large  spontaneous  decay  from  excited  states.  Thus  the  main  heating  is  induced  by  onsite  two-photon 
Raman  transitions  which  do  not  drive  neighboring- site  hopping  but  convert  the  energy  difference  between  two 
Raman  photons  to  mechanical  energy  of  the  lattice  system  [13, 14].  Note  that  in  the  present  optical  Raman 
lattice,  due  to  the  antisymmetry  of  the  Raman  potentials  the  onsite  intraband  scattering  (s  s  bands)  is 

forbidden,  and  only  the  interband  scattering  (s  p  bands)  can  heat  up  the  system.  This  distinguishes 

essentially  from  the  conventional  schemes  which  apply  plane-wave  and  red-detuned  Raman  beams  and  have 
both  inter-  and  intraband  onsite  transitions  [11-15].  The  life  time  of  the  trapped  cold  atoms  can  be  estimated  by 
calculating  the  change  rate  of  the  mean-mechanical  energy  of  an  atom.  For  comparison,  we  consider  both  the 
present  optical  Raman  lattice  system  (with  the  heating  rate  denoted  as  dEOR/ d t)  and  the  conventional  laser- 
assisted  schemes  (denoted  as  dECo/dt).  The  rates  of  change  of  the  mean-mechanical  energy  are  given  by 

^  =  Vr  Y>* (k>  k')^  =  rOR>  (18) 

dt  N  ky 


dEco 

dt 


1 

N 


Y^wss( k,  k')  +  Y^wsp( k,  k') 

k,k'  k,k' 


8lo  —  Tco> 


(19) 


where  N  is  the  number  of  lattice  sites,  and  wss  (wsp)  represents  the  s-s  ( s-p )  band  scattering  rate,  obtained  by  the 
product  of  the  two-photon  Rabi- frequencies  Qss  (Qsp)  and  the  transition  probabilities.  In  the  above  formulae  we 
have  denoted  by  T0r  and  Tco  the  heating  rates  of  the  optical  Raman  lattice  system  and  the  conventional  lattice 
systems,  respectively.  Note  that  the  s-s  band  onsite  transition,  e.g.  from  an  initial  state  with  momentum  k  to  the 
final  state  k'  has  the  two-photon  detuning  As  (k,  k')  =  8lo  +  Ek  s  —  Ek\s,  with  Ek  s  the  s-band  spectrum. 
Similarly,  for  s-p  band  onsite  transition  the  corresponding  two-photon  detuning  reads 

Ap(k,  k')  =  Ek'iP  —  EkyS  —  8cj,  with  Ek\p  the  energy  spectrum  for  the  p-band  states.  Typically  the  bandwidths 
of  s  andp  bands  are  much  less  than  the  s-p  band  gap  which  is  Esp  =  2  ( VoEr)1/2.  We  then  approximate  that 
As(k,  k')  Pd  Scjand  Ap(k,  k')  «  Esp  —  fe,  and 


_ 4  iaPi3 _ 

(2| Esp  -  M  +  2  \nsp\2/\Esp  -  8u |)2  +  4  \Qsp\2  ' 


(20) 


wss 


_ 4  ia5i3 _ 

(2iM  +  2ia5i2/iM)2  +  4ia5i2‘ 


(21) 


Here  Qss  and  Qsp  are  the  two-photon  Rabi- frequencies  for  the  intraband  and  interband  transitions,  respectively. 
For  the  parameter  with  =  5ER,  one  can  verify  that  Qsp  ^  0.5f)ss.  Substituting  this  result  into  the  equation  of 
dEco/dt  we  solve  numerically  that  the  minimum  heating  rate  for  rCo  is  obtained  by  setting  8uo  «  0AEsp. 
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Figure  4.  (a)  Sketch  of  the  spin- 1/2  two  copy  version  of  the  QAH  system.  The  system  is  driven  into  the  h -J2~K  spin  model  with  large 
repulsive  Hubbard-  U.  (b)  The  phase  diagram  versus  J2/ A  and  K/J\.  A  large  region  of  chiral  spin  liquid  phase  is  obtained  due  to  the 
three-spin  interactions.  The  green  point  corresponds  to  7i  =  2 J2  «  3.2 K. 

Namely,  for  the  conventional  laser- assisted-hopping  schemes,  the  minimum  heating  (denoted  as  Tqq)  requires 
that  the  frequency  difference  between  Raman  lasers  be  close  to  the  half  of  the  s-p  band  gap  [13, 14].  With  these 
results  in  mind  we  obtain  directly  from  the  above  two  equations  the  following  relation 


Tor  — 


1  A 

16  Esp  -  A 


rmin 
CO  • 


(22) 


The  life  time  can  be  estimated  by  r  ^  V0/TOr ,  V0/Y{ co-  In  the  present  optical  Raman  lattice  we  can  set  that 
A  C  £sp,  and  then  we  have  T0r  T^o11,  which  shows  that  the  present  optical  Raman  lattice  has  much  less 

heating  than  that  in  the  conventional  schemes.  In  particular,  with  the  parameter  regime  used  in  figure  2,  one 
finds  that  r0R  —  0.03T™  and  for  =  0 .5ER  the  life  time  of  the  optical  Raman  lattice  r  ~  1.67  s,  which  is 
extremely  long  enough  for  realistic  experiments.  We  note  that  other  more  complicated  possible  heating 
mechanisms,  such  as  the  multiband  effects  and  atom- atom  interaction,  are  not  considered  in  the  present  study, 
but  deserves  future  efforts  in  more  systematic  investigations. 


6.  Chiral  spin  liquid  phase 


Finally  we  turn  to  the  realization  of  the  J\-J2-K  model  and  show  it  has  a  large  parameter  region  to  support  the 
highly- sought- after  chiral  spin  liquid  (CSL)  phase  [43-45].  For  this  we  consider  a  spin- 1/2  two  copy  version  of 
the  QAH  model  together  with  repulsive  Fermi  Hubbard  interaction  [figure  4  (a)],  and  have  the  total 
Hamiltonian  that  H  =  Ct(k)7F0(k)C( k)  +  HintwithC(k)  =  (ca  j(k),  4 |(k),  4j(k),  4f(k))r  and  obtain 


Ti  o(k)  =  X  da(k)cra  0  I, 

(23) 

a—x,y,z 

Hint  =  Tij  j , 

(24) 

i 


where  U  is  the  strength  of  Hubbard  interaction.  For  simplicity  we  take  that  mz  =  0.  In  the  single-particle  regime 
each  spin  species  forms  a  QAH  insulator  with  the  same  Chern  number,  while  in  the  large-  U  regime,  the  double 
occupancy  of  each  site  will  be  fully  suppressed,  and  the  system  becomes  a  Mott  insulator.  We  can  then  derive  an 
effective  spin-model  by  considering  the  perturbation  expansion  with  respect  to  t0/U,  t^/U,  with  f0>  t'fl  being 
small  compared  with  U  [46, 47].  For  this  we  consider  the  Hilbert  space  with  single  occupancy  and  treat  the 
hopping  term  Tl 0  as  perturbations.  The  effective  Hamiltonian  is  calculated  by 


({^IHeffK^}}  =X 

n 


({v}\nn0\{a'}) 

Un~i 


(25) 


where  |  { a] )  and  |  { a) )  are  two  different  spin  configurations  with  only  single  occupancy.  Note  that  to  reflect  the 
broken  time-reversal  symmetry  in  the  spin-model,  we  should  at  least  keep  the  terms  up  to  third  order  of 
t0/ 1/,  t'^/U,  which  gives  the  three-spin  interactions  through  triangular  loops.  We  then  reach  the  following 
effective  Hamiltonian  for  the  spin  degree  of  freedom 


Heff  =  X/isi-s;+  E^Si-Sj 

(i>j)  ({«.;» 

+  K  sin  (0^)  S;  •  (S,  x  S*),  (26) 

i,j,ke  A 


where  Ji  =  4 t02/H,  J2  =  Ytl/U  with  tx  =  t'A  «  t^,  K  =  24t02t//f/2,  and  ^  is  the  Aharonov-Bohm  phase 
acquired  by  hopping  through  a  closed  triangular  loop.  The  spin  operators  are  defined  by 

Si,z  =  (CiJfC;,  t  -  ciV<>3/2>  Si,x  =  +  ctici,d/2’  and  Si>y  =  -  c,^  c(i  T)/2 .  It  is  clear  that  the  third 
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order  K-  term  emerges  due  to  the  time- reversal- symmetry  breaking.  The  summation  in  the  third  term  means 
that  each  set  of  (z,  j,  k )  consists  of  a  minimum  triangular.  It  can  be  verified  that  <p-k  =  7r/2  when  <fi0  =  7r/4.  In 
this  case  all  spins  experience  a  uniform  magnetic  field  and  the  spin  system  respects  the  emergent  translational 
symmetry  which,  however,  is  not  respected  by  the  original  free  fermion  system.  We  note  that  the  next-order 
coupling  are  four-spin  interacting  terms,  with  the  four  spins  located  in  the  four  sites  of  a  plaquette.  Since  the  flux 
across  each  plaquette  is  i r,  the  four- spin  interacting  terms  do  not  break  time- reversal  symmetry,  and  are  expected 
to  have  much  weaker  effect  on  the  chiral  spin  liquid  phase.  The  magnitudes  of  Jy2  are  fully  controllable  by 
tuning  t0, i  through  V0  and  Raman  potentials. 

We  solve  the  spinon  mean- field  phase  diagram,  as  shown  in  figure  4(b)  (details  can  be  found  in  the 
appendix).  It  can  be  read  that  three  different  phases  are  clearly  dominated  by  different  interacting  terms  in  He ff . 
The  antiferromagnetic  (Neel)  or  stripe  order  is  obtained  when  the  Jx-  or  J2- term  dominates.  In  the  stripe  phase 
the  staggered  spin  order  exists  only  in  the  x  or  y  direction.  On  the  other  hand,  when  the  three-spin  interactions 
(K- terms)  dominate,  the  CSL  phase  results  [48—  ].  In  this  regime,  no  symmetry-breaking  order  exists  and  the 

spin  degree  of  freedom  is  captured  by  the  bosonic  v  =  1  /2  Laughlin  state  which  has  bulk  semion  excitations 
and  chiral  gapless  spinons  in  the  edge  [52].  When  K  =  0,  the  transition  between  Neel  and  stripe  orders  occurs  at 
J2  =  /i/2,  where  the  system  becomes  most  frustrated.  Then  increasing!^ above  the  point  Ji  =  2 J2  «  3.2K  can 
soon  drive  and  stabilize  the  CSL  phase.  Note  that  the  single-particle  QAHE  in  the  present  system  with  t0  =  V2 1\ 
(corresponding  to  /i  =  2/2  in  the  interacting  regime)  has  a  small  bandwidth  that  is  nearly  1/5  of  the  band  gap. 
This  implies  that  one  can  reach  the  Mott  insulating  regime  with  a  Hubbard  U  which  is  not  too  large  with 
respect  to  hopping  coefficients.  Taking  Vq  =  4\^  =  4 ER,  Vr  =  0.52TRandf/  =  7.6£width,  we  find 
/i  ~  2/2  2i r  x  19  Hz  and  K  «  0.53/i,  which  is  in  the  CSL  phase  region.  We  note  that  the  spinon  mean-field 

calculation  only  shows  qualitative  results  of  the  predicted  phases,  while  recent  numerical  simulation  using 
density  matrix  renormalization  group  method  also  confirms  the  CSL  phase  in  a  similar  spin  model  [  ] .  More 

advanced  investigations  of  the  current  h~J2-K  model  are  necessary  and  will  be  presented  in  the  next  publication. 

7.  Conclusions 

In  conclusion,  we  have  introduced  the  model  of  optical  Raman  lattice  without  SO  coupling  to  observe  chiral 
topological  phases  for  cold  atoms.  We  predict  the  QAH  effect  with  a  large  gap-bandwidth  ratio  in  the  single¬ 
particle  regime,  and  in  the  interacting  regime  we  realize  the }\-}2-K model  which  supports  the  chiral  spin  liquid 
phase.  The  QAH  state  can  be  detected  with  a  minimal  measurement  strategy,  say,  by  only  measuring  the  Bloch 
states  in  the  two  symmetric  momentum  points  of  the  FBZ.  The  minimized  heating  in  our  scheme  and  vast 
tunability  in  parameters  imply  high  feasibility  for  the  observation  of  both  the  single-particle  and  strongly 
correlated  topological  states.  Generalization  of  the  present  optical  Raman  lattice  scheme  to  other  situations,  e.g. 
the  high-orbital  bands,  3D  systems,  and  more  exotic  lattice  configurations,  shall  realize  different  classes  of 
topological  states  which  might  even  have  no  prior  analogue  in  solids.  Especially,  the  correlation  effects  on  such 
topological  phases  should  be  particularly  interesting.  This  work  opens  a  broad  avenue  in  both  theory  and 
experiment  for  the  studies  of  exotic  topological  states  with  cold  atoms. 
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Appendix 

A-l.  Optical  raman  lattice 

As  described  in  figure  1  of  the  main  text,  the  electric  field  of  the  in-plane  incident  laser  beam  is 

E(x)  =  Eo(cosay  +  sin  o:z)el(/cl*_a;t)  +  Eqzc12^-^.  (Al) 

The  initial  relative  phase  between  the  light  components  of  frequencies  u  and  2uj  is  irrelevant  for  the  present 
study  and  is  neglected.  The  in-plane  polarized  components  (x  and  y  polarization  components)  generate  the 
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standing  waves 

Exy  =  2E0cosa[cos(kiX  —  0  —  (j)  +  6(f))y  +  cos(k\y  +  0)x\e^~u;t+e+^~6(^\  (A2) 


Here  6  is  the  phase  acquired  through  the  path  from  mirror  M3  to  lattice  center,  0  represents  the  phase  acquired 
by  the  laser  beam  propagating  along  the  path  from  lattice  center  to  the  mirror  Ml5  then  to  M2,  and  finally  to  the 
lattice  center  again  (refer  to  figure  1  of  the  main  text),  and  6(f)  is  the  phase  tuned  through  the  electric- optic 
modulator.  For  our  purpose  we  set  that  6<f)  =  tt.  On  the  other  hand,  the  out-of- plane  light  components 
generates  the  standing  wave  as 


Ezz  =  4  E0  cos  L(x  +  y)  —  <fi  + 


+4 E0  sin  a  cos  —  (x  +  y)  — 


j  cos  |~ 

4>  -  64> 


h(x  2^  +  y 


A 


i2  —  <jjt+6+4>—  - 


COS 


[l(*  -  '>  - 


0-^0  _  0  _  £ 


1  Jz 

i^-a;f+0+0-£0+y)  ^ 

(A3) 


where  the  additional  it/ 2 -phase  shift  in  the  later  term  is  due  to  the  1  / 4-wave  plate  for  the  z -component  light 
with  frequency  uj  placed  in  the  path  from  mirror  M3  to  lattice  center  (figure  1  (a)  of  the  main  text).  Note  that 
there  is  no  interference  between  the  two  components  with  frequencies  uo  and  2uo.  It  can  be  verified  that  the 
magnitudes  of  the  phases  ( 6 ,  (j)  —  6(f))  only  lead  to  global  shift  of  the  lattice,  and  therefore  are  irrelevant  to  our 
present  study  (they  also  do  not  affect  the  relative  phase  </x  —  (f)y  relating  to  the  Raman  potentials  V#*  and  VRy). 
We  then  set  these  phase  factors  as  zero  to  facilitate  the  description. 

Together  with  the  periodic  Raman  potentials  induced  through  both  the  in-plane  laser  and  the  one 
propagating  in  z  direction  and  having  frequency  uj  —  6uj,  the  total  effective  Hamiltonian  for  the  optical  Raman 
lattice  is  given  by 

p2  -f  p 2 

H  =  — - E  +  ^sq(x,  y)  +  V^ix,  t)  +  VRy(y,  t),  (A4) 

2m 


where  m  is  atom  mass,  the  double- well  square  lattice  potential  V^q  (x,  y),  the  time- dependent  Raman  potentials 
VRx  ( x ,  t)  and  VRy  (y,  t)  take  the  forms 


Ki(x,  y)  =  M)[cos2(fc!x)  +  cos2(fciy)]  +  V0sm2[k(x  +  y) ]sin2[/q(x  -  y )] 


+Visin2  /(x  +  y) 


1  2\h 

cosz  — 

J  12 


—  (x  -  y) 


\4c(x,  t)  =  2V#  cos (6u)t  —  (j) y)cos(/cix),  VRy(y ,  t)  =  2V#cos (6u;t  —  (f)x)cos(kiy ). 


(A5) 

(A6) 


The  third  term  in  V^q  (x,  y)  leads  to  an  onsite  energy  offset  A  =  V\  between  A  and  B  sites,  with  V\  being  small 
compared  with  Vq.  The  second  term  with  amplitude  reduces  the  difference  in  height  of  the  barriers  along  the 

AB- bond  and  the  diagonal  ( AA/BB )  directions.  Thus  it  can  enhance  the  diagonal  tunneling  relative  to  the 
hopping  coupling  between  A  and  B  sites,  providing  vast  tunability  in  parameters.  The  neighboring  hoppings 
between  A  and  B  sites  are  restored  when  6u  «  A.  Here  we  consider  only  the  s-orbitals  ^^(r)  (p  =  A,  B ), 
which  are  of  even  parity. 


A-2.  Spinon  mean  field  approach 

We  solve  the  spinon  mean-field  phase  diagram  for  the  J1-J2-K model  with  (f>-k  =  it  j2.  This  model  contains  at 
least  three  phases.  First,  when  Jx  dominates,  the  system  is  unfrustrated  and  it  supports  a  Neel  anti- ferromagnetic 
order.  Secondly,  when  J2  dominates,  the  system  is  also  unfrustrated  and  can  have  a  stripe  anti- ferromagnetic 
order,  in  which  case  the  staggered  spin  order  exists  only  in  the  x  or  y  direction.  Finally,  when  K  is  large  enough, 
the  system  prefers  a  chiral  spin  liquid  state  [43, 44]. 

The  different  phases  can  be  studied  using  trial  (mean  field)  wave  function  method.  We  introduce  the  anyonic 
spinons  f.  =  f^)T  to  represent  the  spin  operators  as  St-  =  f?  —f{  under  the  particle  number  constraint 

=  1.  In  the  mean  field  theory  for  £7(1)  spin  liquid,  the  spin  interactions  can  be  rewritten  as  the  following, 


5,  •  S, 


"  2  AijXji’ 


Si  x  Sj  •  Sk  — 


2  i  6 


{[XijXjkXki  +  cyclic (z)7c)]  -  h.c.}, 


(A7) 

(A8) 


where  x~  —  xjj  =  is  the  spinon  hopping  operator.  Here  we  do  not  consider  the  spinon  pairings  since  our 

interest  is  mainly  focused  on  the  £7(1)  chiral  spin  liquid  phases.  In  general  the  spinon  hopping  term  is  complex, 
and  the  spin  chirality  term  can  give  rise  to  a  phase  e1(^A,  with  (f)A  =  Arg((x;z)  (Xkj)  {Xik) )  the  ^ux  experience  by 
spinons  after  hopping  through  a  close  triangular  loop.  For  chiral  spin  liquid  state,  in  this  work  we  have 
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numerically  verified  that  the  ground  state  corresponds  to  <pA  =  y .  Actually,  with  a  Landau  gauge  choice  for  the 
mean  field  theory  one  can  confirm  that  the  i r / 2  -  flux  state  in  each  triangle  has  the  lowest  energy.  Therefore,  for 
convenience  we  introduce  the  trial  mean  field  parameters 

Xl  =  (Xn+l*)*  Uj=odd  =  —  (X«+lx)  Uj=even 
(X«+ly)  l>;=odd  —  (X«+ly)  ly—even* 

X2  =  (Xii+1,+1,)  =  X* 

to  decouple  the  spin  interactions.  Here  x2  is  assumed  to  be  real  and  the  hopping  phase  is  carried  by  X\  f°r  the 
chiral  spin  liquid  phase.  With  other  gauge  choice  for  the  mean  field  parameters  we  shall  get  the  same  phase 
diagram.  It  is  easiy  to  see  that  with  the  above  trial  mean  field  parameters  the  spinons  experience  a  uniform  1/(1) 
gauge  field,  with  the  magnetic  flux  through  each  triangular  being  7r/2  and  through  each  plaquette  being  7r[44] . 

Furthermore,  since  the  system  may  contains  symmetry  breaking  orders  at  some  parameter  region,  we  also 
need  to  introduce  the  following  magnetization  order  parameters  to  describe  the  Neel  order  and  stripe  order, 
respectively 

Mn  =  (-1  )*■■+* (Sf),  Ms  =  (-l)MSf).  (A9) 


Now  we  can  decouple  the  spin  interactions  by  these  (trial)  mean  field  parameters  as 


Si  •  S*+lJjCi=odd  —  ^ 

‘  Si+lJ*i=even  —  I 


(-“XiXii+1*  +  h.cj  +  i  |Xl|2  +  ( —  1  )Xi+yi (M„ Sf+  ix  -  MnS?)  +  M2, 

(-^xtXu+ix  +  h-cj  +  1  |Xl|2  +  (-1  )x‘+y-(MnSf+lx  -  MnS?)  +  Ml 

Si  •  Si+ly\yi=odd  =  (-^XfX,I+ly  +  h.cj  +  i  |Xl|2  +  (-1  y*X(MaS?+ix  -  MnS?)  +  Ml 

S,  ■  Si+ly |yi=even  =  -  (-iXlXii+lr  +  h.cj  +  i  |Xl|2  +  (-l)^(M„Sf+li  -  M„Sf)  +  M2, 

S,  •  Si+Il+Iy  =  (-tX2X,+  li+ly  +  h.cj  +  i  |X2p  +  (-1  r-+y.(MnS?+lx+y  +  MnS?)  -  M2n 


+  (— 1  )x>(MsS?+lx+ly  -  MsS ?)  +  Ml 

Si  x  Sj  •  Sk  =  ^-[(XijXjk)Xki  +  cyclic ( ijk)  -  h.c.]  -  2-^[|X2X2|e-^A  -  h.c.] 
4 1  J  J  4 1 


1 

4  i 


lxfx2le 

(Xki) 


Xki  +  cyclic  (ijk)  -  h.c. 


-  2-t[|Xl2X2|e  '■&  -  h.c.]. 
4i 


Notice  that  both  the  Neel  order  Mn  and  the  stripe  order  Ms  are  collinear,  so  they  don’t  appear  in  decoupling  the 
spin  chirality  interaction  term  Si  x  Sj  •  Sk.  In  the  chiral  spin  liquid  phase  Mn  =  Ms  =  0,  and  from  the  above 
expressions  we  find  that  the  spinons  experience  a  uniform  magnetic  field  which  leads  to  the  quantum  Hall  effect 
for  spin  degree  of  freedom.  Then  the  ground  state  of  the  chiral  spin  liquid  phase  is  captured  by  the  bosonic 
v  —  1/2  Laughlin  state  [52]  which  has  chiral  gapless  anyonic  spinon  excitations  in  the  edge  [43, 44].  The  phase 
diagram  can  solved  self- consistently  based  on  the  above  decoupled  formulas. 
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